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Heavy ion collisions at RHIC/LHC energies are well described by the (nearly ideal) hydrodynam- 
ics. Last year this success has been extended to higher angular harmonics, v n ,n — 3.. 9 induced by 
initial-state perturbations, in analogy to cosmic microwave background fluctuations. Here we use 
hydrodynamics to study sound propagation emitted by quenched jets. We use the so called "geo- 
metric acoustics" to follow the sound propagation , on top of the expanding fireball. The conical 
waves, known as "Mach cones", turn out to be strongly distorted. We show that large radial flow 
makes the observed particle spectra to be determined mostlly by the vicinity of their intersection 
with the fireball's space-like and time-like freezeout surfaces. We further show how the waves modify 
the freezeout surfaces and spectra. We end up comparing our calculations to the two-particle cor- 
relation functions at RHIC, while emphasizing that studies of dijet events observed at LHC should 
provide much better test of our theory. 

PACS numbers: 



I. INTRODUCTION 

A. Hydrodynamics and sounds 

Jet-induced correlations have short and remarkable 
history. Two-hadron correlation functions measured in 
AuAu collisions at RHIC have revealed structures known 
as "ridge" and "shoulders" (about 2 rad away from trig- 
ger) These structures were originally belived to be as- 
sociated with jets, but have recently been explained as 
hydrodynamical "harmonic flows" , basically a sound cir- 
cles created by initial local perturbations. This sub- 
ject has been discussed extensively in three papers [T]- 
3 coauthored by one of us, which have very extensive 
introduction and references, to which the reader may re- 
fer. The first paper of this series [TJ has been discussing 
some generic features of the sound propagation and 
possible role of electric (dual-magnetohydrodynamical) 
corona. The second paper [2] has included more detailed 
studies of the phases of the perturbations and possi- 
ble ways of their experimental/theoretical studies. The 
third paper [3j used the so called Gubser flow, which al- 
lows analytic treatment of all angular harmonics, with 
and without viscosity. (One element of this last paper, 
the modification of the freezeout surface by the sound 
waves, will be extensively used at the end of this work.) 
While thermo/hydrodynamics was traditionally expected 
to work for majority of soft particles only, it describes 
well these harmonic flows including the so called "inter- 
mediate transverse momenta" region, p t = 2 — 3 GeV . 
Such particles are rare, about one such particle per 
event among thousands. So, once again, one finds that 
thermo/hydrodynamics is very robust and can describe 
not only the behavior of an average particle, but also 
the shape of rather far-reaching tails of thermal/hydro 
spectra. 

These developments had (perhaps temporarily) take 
attention away from the issue of jets and their interac- 



tion with the medium. However, as soon as the trig- 
ger hadron has pt above the intermediate pt domain, 
say pt > 6 GeV, it gets out of reach for collective flows 
and comes predominantly from hard collisions, or "jets" . 
So, the experimental input to the present paper are two- 
particle correlations in which the trigger is "hard" , while 
one (or more) associate particles are in the hydro do- 
main, p t = 1 — 3 GeV (hard-intermediate correlations, 
for short). 

The idea, that once the energy is deposited into the 
medium by a jet will be resulting in sound perturbations 
in the shape of the Mach cone, has been proposed in Refs 
[HE]. (It may resemble similar idea discussed in 1970's 
for non-relativistic nuclear collisions. However it did not 
work, because nuclear matter is dilute and not a partic- 
ularly good liquid, unlike sQGP under discussion now.) 
In the context of strongly coupled QGP, the problem has 
been addressed in the framework of AdS/CFT. As de- 
tailed by Chesler and Yaffe [5], and Gubser et al [7], the 
stress tensor solution obtained by holographic imaging is 
in remarkably good agreement with hydrodynamical so- 
lution detailed in [4], displaying the Mach cone in all its 
glory. Since the time-averaged sound velocity over the 
QGP, mixed and hadronic phases is < c s > s» 0.4 , the 
expected Mach cone angle 

9 M = arccos [ < ° s > ) w 1.159 rad w 66.4° (1.1) 

V "jet / 

roughly matches the angular positions of the "shoulders" 
at (f> sa 7T ± 9m in the correlation functions. However, as 
the reader will see from what follows, a flow normal to 
the Mach cone is only one of several geometrical issues 
involved. Another, equally important effect, is the radial 
flow of extra matter produced by the wave. As we will 
see, it is defined mostly at the point where the sound 
cones intersect with the time-like and space-like freezeout 
surfaces. 

One more geometrical effects come from the interplay 
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between the jet "stopping distance" and the actual size 
of the fireball (along the jet direction). Obviously, when 
the two are equal, the amplitude of the wave is at its 
maximum. This effect has been discussed before in direct 
hydrodynamical simulations of the cones, such as, e.g., 
that done by Betz et al [8], which we, to some extend, 
will follow. 

The main technical difficulty of the problem is in cor- 
rect treatment of the sound propagation on top of the 
expanding fireball. This problem has been first discussed 
by Casalderrey-Solana and Shuryak in [24] which mod- 
eled the fireball expansion by a Hubble- like overall expan- 
sion of the metric, the same Friedman-Robertson- Walker 
metric as used for Big Bang. (Of course, it creates new 
space rather than expanding, so it is not literally corre- 
sponding to the problem at hand, it was just a technical 
first step.) The physics focus of that paper was the effect 
of time-dependent sound velocity, especially if the phase 
transition is of the 1st order and can vanish at some in- 
terval of temperature T . The interesting finding was a 
creation of secondary - and convergent - sound waves. 
This idea was further discussed in ;1 in connection with 
the "soft ridge" issue. However, if the current lattice data 
on the speed of sound are correct, the calculated effect of 
the reflected wave was shown to be too small to explain 
the "hard ridge" . 

In this paper we developed a different method based 
on "geometric acoustics. 

Let us now explain the structure of the paper. In Sec 
II we work out the equations of the propagating sound 
in a moving fluid using the geometric acoustics [5]- We 
consider both relativistic and non-relativistic cases. Then 
in Sec III we construct the Mach cones fomed by sound 
waves emitted by the associated jet, and circles fomed 
by sound waves emitted by a jet which is originated at 
the same hard collision point but is different from the 
trigger and associated jets. Meanwhile, all the cones and 
circles are solutions of the relativistic and non-relativistic 
equations. In Sec IV we proceed to discussion of the issue 
related to intersection of the Mach cone surface with the 
fireball's freezeout surface in 3D. Finally, in Sec V we 
calculate the spectra of secondaries comparing them with 
the experiment as well. We will use our assumption that 
most of the contribution to the spectra come from the 
region which is the "edge" formed by the intersection of 
the surfaces. 



B. Comments on the jet quenching 

This complicated phenomenon, suggested by Bjorken 
in 1982, has very long history, which we of course would 
not go into here. For the purpose of this work it is 
not important to know the microscopic mechanism of jet 
quenching, and for modelling we will use simple forms of 
the energy dE/dx deposited into ambient matter. Thus 
this subsection will only include comments relevant for 
what follows below. 



Hard collision events of partons are a particular type 
of "initial state perturbation". It is important [10] that 
those generates four (not two!) jets. Two of them ap- 
proximately balance their large transverse momenta One 
of these two jets is the parent of hard "trigger" hadron, 
the other called the associate jet, propagates through the 
medium and deposits certain amount of energy into the 
ambient matter. Below we will follow what happens with 
it, provided it can be described hydrodynamically. 

Two more bremstrahlung cones, or jets, are directed 
forward and backward in the beam direction. Since these 
ones do not have large transverse momenta, they are 
rarely discussed, if at all. However, they do produce ex- 
tra multiplicity, roughly the same as two other jets, which 
should be seen in the correlation functions. As suggested 
in [TO], those would naturally explain the "hard ridge", 
an extra structure near the triger direction in the az- 
imuthal angle, but long-range in the rapidity. In [1] the 
main puzzle has been discussed: why this perturbation 
remains correlated with the trigger direction rather than 
create a sound circle and move ±1 radian away, creating 
"shoulders" instead? The suggested answer was flux tube 
formation, natural if hydrodynamics is elevated to dual 
magnetohydrodynamics. In this paper we will however 
not discuss that option, restraining to the usual sound 
waves. 

For modelling energy loss one needs to know how the 
jet energy loss depends on the matter (entropy) density 
s(x) as well as the length L travelled by the jet in the 
matter from the moment of its production. 

The simplest assumptions, resembling what happens 
in QED, is that it is proportional to matter density and 
is independent of x, —dE/dx ~ s * const(x). It did not 
however took long to notice [HJ [12] that this assumption 
is in conflict with the experimental data, on jet quench- 
ing dependence on the azimuthal angle 4> relative to the 
impact parameter direction. Those ususally are charac- 
terized by the parameter 

v 2 (pt) =<cos(20) >, (1.2) 

(not to be confused with the elliptic flow: we now hard 
transverse momenta p t > 6GeV). This quantity can 
also be seen as related to the difference between the 
hard hadron spectra from jets travelling in the x direc- 
tion (<f> — 0, in-the-reaction-plane) and in the y direc- 
tion ((f) — 90°, out-of-plain) . Not going into the his- 
tory of the V2 debate, let us comment on its current sta- 
tus. There are two ideas on the market, reproducing the 
vi(vt > GGeV) data: 

(i) One [T3] is that the quenching is not proportional to 
density but is enhanced in the near-T c region or the so 
called "mixed phase" . 

(ii) Another is based on the strong coupling theory 
(AdS/CFT) picture of matter, resulting in the energy 

loss on Eg 

~. | strong coupling ^ T X , (1-^) 
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While this regime may look similar to the weak coupling 
BDMPS result 



clE 
dx 



weak coupling ^ T X 



(1.4) 



it is in fact qualitatively different. To see this one should 
recall that T ~ 1/x 1 ^ 3 due to longitudinal expansion. 
As a result, the strong coupling result diverges at large 
x, while the weak coupling one is spread over the whole 
path. Thus the former one does describe the v 2 data, 
while the latter one does not. 

Two proposed explanations are very different in na- 
ture, the near-T c "blackness" of matter [T3] versus phys- 
ical growth of the jet "falling into the IR" in O ITS] . 
Yet it is hard to tell them apart in practice, because at 
RHIC the near-T c region corresponds to the proper time 
t = 5 — 9 fm/c which is very close to the typical time it 
takes for the associate jet to traverse the fireball. Perhaps 
one would need to do a detailed comparison of RHIC and 
LHC results to see which one describes the data better, 
as at LHC the near-T c region is shifted to later times. 

Another feature of the strong coupling results is that 
large fraction of energy is released at the stopping point, 
see details in [T5]. In our calculations below, there is one 
more reason for the importance of the endpoint of the 
jet path inside the matter: the less time perturbation in 
matter has to travel, the larger sound amplitude it has. 



II. SOUNDS IN A MOVING FLUID AND THE 
GEOMETRIC ACOUSTICS 

As detailed in hydrodynamics textbooks (e.g., by 
Landau-Lifshitz [5]) the rays describing the sound prop- 
agation can be described in the "geometric acoustics" 
approximation which uses the analogy between the 
Hamilton- Jacobi equation for the particle and the wave 
sound equation. The velocity and direction of the 
sound propagation can be derived with the help of the 
eikonal function. The Hamilton equations of motion for 
"phonons" (or "particles" of the sound) have the follow- 
ing generic form 



df du{k 1 r) 
dt ~ dk ' 
dk duj(k,f) 
dt df 



(2.1) 
(2.2) 



driven by the (position dependent) dispersion relation 
uj(k,f). For clarity, let us start with the simplest non- 
relativistic case, namely of small velocity of the flow, u <C 
1. In this case the dispersion relation is obtained from 
that in the fluid at rest by a Galilean transformation, so 
that 



w(k, r) = c s k + ku . 



(2.3) 



In the trivial case of the constant flow vector, u — 
const(r), the first of the above equations just obtains 



an additive correction 
df 



-=c s n % + u, 



(2.4) 



where ftj: = k/k is the unit vector in the direction of the 
momentum. The second equation remains trivial 



dk 



0. 



(2.5) 



as there is no coordinate dependence anywhere. So, a 
point perturbation in a moving fluid simply produces a 
moving sound circle, motion with the speed of sound, 
then the position 

After this "warm-up", let us do a non-trivial problem, 
simple enough to be analytically solved. Let us consider 
a (generalized) Hubble-likc flow 



Ui{r) = HijTj 



(2.6) 



with some time and coordinate independent Hubble ten- 



sor Hij. The equation (2.2) now reads 



dki 
~dt 



— kj 



(2.7) 



So, if the Hubble tensor is symmetric and can be diag- 
onalized with real eigenvalues Hi, H 2 , H3, then the so- 
lution in its eigenframe is exponential momentum con- 
traction ki(t) = exp(—Hit)ki(Q). However, if the Hubble 
tensor contains an anti-symmetric part, the eigenvalues 
can be imaginary, so that the k is rotating around the 
vector CijkHjk. 



The first equation (2.1) reads : 



dn 
lit 



c s nAt) + Hijrj (t) 



(2. 



In the simplest case when the Hubble matrix is propor- 
tional to the unit matrix the solution is simply a linear 
addition of the Hubble expansion and sound motion: 



fit) 



c s tn^ 



f(0) exp(+Jft) 



(2.9) 



This example can be applied to real hydro explosion, in 
which the late-time flow can indeed be well approximated 
by a Hubble form. For central collisions we will only 
discuss, the Hubble tensor is even isotropic. The only 
complications is that this pattern needs some time to be 
developed, so below we will use it only later than some 
time to m 5 fm [T7] . So the two-dimensional form of the 



equations (2.1) and (2.2) to be used is 



dk\ 
lit 
dk 2 
lit 



dri 

~dt 
dr 2 
~dt 



-Hki(t)Q(t-t ) 

-Hk 2 (t)e{t-t Q ) 

c s cos (0) + Hn(t) 9(t - t ) 
c s sm(6)+Hr 2 {t)e{t-t ). (2.10) 
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The initial conditions to this system are specified as fol- 
lows: 



fci(O) = cos (0), fe 2 (0) =sin(0) 
n(0) = x(0), r 2 (0)=y(0), 



(2-11) 



where x(0) and y(0) are constants and measured by 
fm 's. They are the coordinates of the hard collision 
point where a jet is produced and simultaneously emits 
phonons. For the angle 9 we take the value of the Mach 
angle relative to the jet velocity from equation (1.1): 
cos (9) = cos (9 M ) =< c s > /v jet . 

The corresponding analytical solutions will be 



h(t) 
k 2 (t) 
ri(t) 



= cos(9)e- m@ ^-^ 
= sm(9)e- m@ ^-^ 
= [c s cos (9) t + x(0)] Q(t - t ) + 

c s cos (9) (1 - e^*"* ) - 5H e ff (*-*°)) 



H 



ra(t) 



-x{0)e H(t - to > (@(t-t )-l) 

c s sm(9)t + y(0)]&(t-t ) + 
c s sin (9) (1 - e"^-^ - 5H e H{ - t - f ^) 
H 

■y(0)e H ^] (e(t-to)-l). (2-12) 



A bunch of rays emitted to all directions from some 
initial point will thus be at a surface of an expanding 
sphere, as in constant flow: 



(r(i)-f(O) eM+Ht)f 



(2-13) 



with a sphere center moving "with a flow" . 



In heavy ion collisions the explosion is relativistic, with 
transverse velocity reaching u ~ 0.7c. Therefore the lo- 
cal Galilean transformation should be substituted by the 
Lorentz one, which gives 



to(k, r) — c s k cosh(Y) + ku sinh(Y) 



(2-14) 



where u means a unit vector in the direction of u, and Y 
is the rapidity. Now it is no longer possible to separate 



the equations (2.1 1 and (2.2 1, and they should be solved 



together numerica ly. Su 



DStitution of the equation (2.14) 



into (2.1) and (2.2) gives us the following system of four 



differential (relativistic) equations: 




FIG. 1: The a is the angle between the associated jet's mo- 
mentum and the axis x. The arrow inside the circle is the 
momentum of the associated jet which originates from the 
hard collision point opposite to the trigger jet (not shown) in 
the x direction. 



dki 

lit 




ki 



dn \ ^l-m(ri + ri)e(t-t , ) 

d ( H ri e(t-t ) 

d ( Hr 2 Q(t-t ) 
On [y/l-HZtf + rfMt-to) 



dk 2 




ki 



dri 

~dt 



dr 2 \ y/l-H*(r? + r*)Q(t-t 0) , 

d ( Hrx®(t-to) 
dr 2 { y/l-H^rf + rZMt-to), 

d ( Hr 2 e(t-t Q ) 
dr 2 [ y /l-H^ + ri)Q(t-t ) i 

fei 1 



yfc?TfcJ v/l - H\rt + r*)Q(t - t ) 

HnQ(t-t ) 



dr 2 
~dt 



y/l-H*(r?+r2)Q(t-t ) 



k- 2 



^/kJTkl v/l - H 2 (r? + r|)0(t - t ) 
Hr 2 Q(t-t ) 



^l-H2(r?+rl)O(t-t ) 



(2.15) 
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FIG. 2: (Color online) A Mach cone originated from a jet moving to the left from the point x(0) = 5.0 fm and 7/(0) = fm 
at the angle a = 0. The path of the jet corresponds to L stop p ing = 6 fm, the radius of the inner circle - R in = 6 fm, and 
the radius of the outer circle, containing the fireball at its highest extension Rj — exp(H-(tf — to)) — 10.14 fm. The cone 
created by the open squares is the soluti on of t he non-relativistic equations (2.101, and the cone created by the open circles is 
the solution of the relativistic equations (2.151. The filled circle shows the hard collision point, and the open diamond shows 
the stopping point of the jet. The values of the angle a shown in each figure, is defined in Fig.l. 



with the same initial conditions as in (2.11 ) 



III. DISTORTED CONES AND CIRCLES OF 
THE EXPANDING FIREBALL 



We solve the systems of equations (2.10) and (2.15) 
at freezeout time tj = 12 fm, speed of sound c s = 0.4, 
and the Hubble constant H = 0.075 fm . The tf is the 
time at which the fireball reaches the freeze-out tempera- 
ture. (In Appendix A, we clarify a value for the freezout 
time used, along with establishing the kinetic freezeout 
temperature.) 



In order to simplify geometry of the collision, we con- 
sider only central Au — Au collisions. The notations are 
explained in Fig.[TJ jets which have different impact pa- 
rameters in respect to the fireball are defined via the 
angle a, between the direction of the associate jets' mo- 
mentum and the axis x in which the trigger jets go. (We 
do not show the trigger jet in Fig.[lJ as a reminder that 
we do not include the sounds from it in what follows, 



G 
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FIG. 3: The Mach cone (open squares) as in Fig.[2] and the ring (open circles) which is originated from the hard collision point 
x(0) — 5.0 fm and y(0) = fm at the angle a = 0. These two structures are the solutions of the non-relativistic equations 
(2T0I. 



considering its energy loss to be negligible. Due to trig- 
ger bias effect, we place hard scattering events at fixed 
distance, 1 fm, from the nuclear edge.) 

Quenching of the associate jet is the source of the 
phonons, emitted from its path along the Mach direc- 
tion. After that they move along the phonon rays as 
detailed in the previous section. We show only solutions 
in the upper part of the transverse plane of the fireball, 
binning angle a into six bins a = 0, 15°, 30°, 45°, 60° 
and 75°, which are depicted in Fig.[|] The case of a = 90° 
does not give rise to any conical structure since the jet 
at this point only strikes the edge of the fireball/inner 
sphere leaving it totally, without making phonons prop- 
agate in the medium. The solutions in the lower part 
are symmetrical to those in the upper part, that's why 



they are not shown. The conical and ring structures are 
the solutions of the systems of equations: non-relativistc 
(I2.10D and relativistic (12 151). 



Note that in all these six figures the jets move to the 
left with the speed of light, from points of their produc- 
tion, located inside the inner circle at a distance 1 fm 
from that circle which in fact is the surface of a nucleus. 
This inner circle corresponds to the case of tf = fm, 
when the firebal/inner sphere initially is not expanded 
as a whole. The outer circle corresponds to the radius 
tf = 12 fm, it shows the maximal extension of the fire- 
ball reached by this time. 

The jet quenching discussed in those figures assume 
the stopping distance of the jets be L stopping = 6 fm, 
shown by the diamonds. ( For comparison purposes 
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we have also calculated cones for the stopping distance 
L stopping = 12 fm not shown, with similar results to 
those in Figs. [2] except for different "rounding" of it does 
not matter much.) This means that we assume the jet 
energy and the energy loss are related by E/\dE/dx\ = 
6/m. The actual values of E,dE/dx are not important 
for these plots and will matter only later, when we will 
determine the actual spectra modification. 

We see that at the angles increasing towards higher 
a, the conical structure created by the phonons has a 
distorted/crooked shape, which is due to the Hubble ex- 
pansion starting from to ~ 5 fm. The jets, produced 
at low and intermediate angles, do not reach the edge 
of the outer sphere, meaning that they do not leave the 
fireball. However, at the higher angles, some jets go out 
from the fireball at some values of the time within the 
range of tf = 5 — 12 fm. One conclusion is that in the 
fireball with the given kincmatical and dynamical con- 
ditions, the shapes of the Mach cones do not differ from 
each other significantly obtained from the non-relativistic 
equations ( |2. 10 1, and from relativistic equations (2.15). 
Therefore, hereinafter we can consider only one of these 
cases, namely one can take the non-relativistic case for 
earring out the other calculations. 

We already mentioned above two more jets which are 
also created by a hard collision which move perpendic- 
ularly to the transverse plane along the beam direc- 
tions. These two longitudinal jets are approximately 
rapidity-independent source which also may emit the 
sound waves. Their propagation we calculated using the 
same equations (2.10) but with the initial conditions in 
(2.11), where ranges from to 360°. Such (rapidity- 
independent) sound circles are depicted in Figs. [3] to- 
gether with the cones (which are for one rapidity only, 
namely that of the associate jet). We will not use those 
circles below: in fact we have calculated them in order 
to compare to what is calculated in [5] by a different 
method. 



IV. INTERSECTION OF THE MACH SURFACE 
WITH TIMELIKE AND SPACELIKE 
FREEZEOUT SURFACES 

Standard expression for a particle spectrum, known as 
Cooper-Frye formula |31) , is a thermal spectrum boosted 
by the 4-vector of the flow velocity m m and integrated 
over the 3 — dimensional freezeout surface. 



E 



^dN 
d 3 p 



P_Uv 

T 



(4.1) 



where at and a s are time-like and space-like parts of 
the freezeout surface, respectively. The function / is 
Bose/Fcrmi/Boltzmann distribution, for pions it is 




P u v 
T 



1 



exp(-p ly M l/ /T) ± 1 



(4.2) 



FIG. 4: The schematic shape of the Mach surface in the 
transverse (x, y) plane at z = and fixed time (upper plot), 
as well as its shape in 3d including the proper longitudinal 
time (lower plot). The Mach surface ctm is made of two parts, 
OCAA'T and OCBB'T. This figure is from [30]. 



simply corresponds to the thermal distribution of mat- 
ter inside the fluid cells. (If needed, the function can be 
appended by an anisotropic part corresponding to viscos- 

ity.) 

We will be interested in a special case of this distribu- 
tion in the region where, at one hand, the hydrodynam- 
ics works well, and on the other, the ratio (pt/Tf is as 
large as possible. (As we are only interested in the tail 
of the function /, so that the Boltzmann approximation 
will always be enough.) This ratio, if large, strongly en- 
hances the effect of the small perturbation, the sounds, 
on the particle spectrum. In fact, for p t = 1 GeV and 
Pt = 2 GeV we will be using below, this large ratio time 
small amplitudes of the sound waves is able to produce 
effects of the order 0(1) in the exponent. This interplay 
of large Pt/Tf and small sound wave is one of the core 
ideas of this paper. 

Before we go to details about realistic freezeout sur- 
faces, in the next section and Appendix A, we would like 
to introduce another qualitative idea, this time related 
with the radial flow. The concept has been introduced 
in [30] and is explained in schematic way in Fig.|4j (Note 
that in Fig. [4] it is assumed that the associated jet prop- 
agates to the right, unlike in the conventions above.) It 
shows the picture in 3d, at fixed longitudinal coordinate 
z = 0. The proper time r runs upward, and the upper 
circle schematically represents the time-like part of the 
freezeout surface, at , approximated by the constant time 
surface, t — tf. The lower circle is the so called "initia- 
tion time surface" , and the conical surface connecting the 
two ellipses is our approximation to the space-like part 
of the freezeout surface, a s . The points A, B, A' and 



B' are the intersection of the sound waves with cr 4 , and 
the region between them contains matter affected by the 
waves: outside it is not affected, by hydro causality. The 
points T, C are intersections of the trigger and the com- 
panion jets with a s . The main idea is that the vicinity of 
four points A, B, A' and B' should dominate the spectra 
(modified by the sounds) because the large upper circle is 
where the radial flow is at its maximum. (We will return 
to this idea at the qualitative level in the next section.) 

Let us now promote this figure from 3-dimension to 
all 4, by adding the longitudinal z direction. We cannot 
plot it in full, and therefore limit ourselves to the t = 
tf time slice. Since we have the intercept of 3 surfaces 
(Mach,er t , a s ) in 4d, it must be a 1-dimensional curve. 
It is in fact two ellipses, one having points A, B on it, 
and the other having A', B' . Let us call these two edges 
ec and e T , for the companion (associated) jet and the 
trigger jet, respectively. (Since the trigger-bias forces the 
associated jet to deposit much larger amount of energy, 
the former one has much larger chance to become visible.) 

Let us calculate the extension of these elliptic curves 
into the longitudinal coordinate z direction, or rather its 
proxi, the pseudo-rapidity rj. In Figs. [5] |6j [7] and [8] we 
show correspondingly the edges ec (ellipses) created by 
phonons at tf = 12 fm, emitted by a jet whose momen- 
tum is directed as in Figs. [2] However, now the increase 
of the angle a from the previous jet to the next one is 
Aa = 5°, rather than 15°. Again the case of the jet pro- 
duced at a = 90° is ignored. Besides, we also consider 
phonons (from all the jets under consideration) moving 
in the lower hemisphere of the fireball. The total number 
of all the considered jets (events) is thirty. 
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FIG. 5: (Color online) Two elliptic structures correspond to 
the edges ec produced by phonons at intersection of all the 
three surfaces. These phonons in turn are emitted by an as- 
sociated jet which moves in the x direction (or parallel to it) 
emitting the phonons in x, y and 77 directions. The largest 
ellipse corresponds to the jet with a = 15°, and the smallest 
ellpise corresponds to the jet with a — 85°. 



n 




FIG. 6: (Color online) The same figure as in Fig.|5]but which 
is complete including the contributions of all the jets with 
a = 15° -85°. 
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FIG. 7: (Color online) The same figure as in Fig.|6]but in the 
(?7, y) plane. The largest ellipse corresponds to the jet with 
a = 15°, and the smallest ellpise corresponds to the jet with 
a = 85°. 



V. APPLICATION TO PARTICLE SPECTRA 
WITH THE MODIFIED FREEZEOUT SURFACE 

The temperature and velocity in Cooper-Fry formula 
have complicated space-time dependence derived from 
hydrodynamics. As in [3] we assume that the freezeout 
surface is an isotherm of the form of T(i, r) = Tt. Ac- 
cordingly, the time-like part of the surface, for example, 
can be written as 



^ = (tf(x,y),x,y,r)) . 



(5.1) 



So at the time tf(x,y) the fireb all reaches the freezeout 
temperature Tf. Note that (4.1 1 contains the vector nor- 
mal to the surface. As to the element of the surface, it is 
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FIG. 8: (Color online) The same figure as in Fig. [6] but in the 
(j), x) plane. The largest vertical corresponds to the jet with 
a = 15°, and the smallest vertical corresponds to the jet with 
a = 85°. 



represented as 

dX p = ^/-ge^xp— —dxdydr] 

ii.r (li< (is, 

= —tf dxdydr] + tfdtfdydrj + tfdtfdxdrj , (5.2) 

where g is the determinant of the metric and e^ v \ p is the 
Levi-Civita symbol. 

For p^dY*^ and p p u^ one may use the formulae from 
Ref. [35] which are the following: 



tf(r) rdrd(f>drj,(5.3) 



p^dT, 11 = [m t cosh(ri - y) - 

dt f (r) 



Pt cos( 



dr 



p^ = (m t cosh(?7 - y)u T - p t cos((f> - (t> p )u r ) (5.4) 



where u T = 7 — (1 — u(r) 2 ) x l 2 , u r = u(r)u T (withu(r) 
given by (2.6)), mt — \Jp\ + m 2 (the m n is the pion 
mass). We extract the function tt(r) from hydro solu- 
tions discussed in Appendix A. 

For each p t and the direction of the particle <j) p there 
is a point on the freezeout surface which maximally con- 
tribute to the spectrum. Indeed, the radial flow grows 
with r and enhance the spectrum, but the surface itself 
makes a turn and ends. One can determine the spec- 
trum from the Cooper-Fry formula and compare it to 
e.g. Gaussian approximation around this point. Those 
two are shown in Fig. [9j as a particular example. The 
lesson is that there exist a sharp peak near the fireball's 
edge, at which the quest for the strongest radial flow and 
for the largest amount of surface area (and thus multi- 
plicity) are in the best compromise. We have checked 



that those peaks indeed produce the spectra in agree- 
ment with what hydro papers reported, and with what 
are observed experimentally. This confirms that the tails 
of the spectra we are interested are dominated by the 
vicinity of such points. 




FIG. 9: (Color online) The contribution to the pion spectrum 
(arbitrary units) vs the radial coordinate r (fm). The solid 
(black) curve is the Cooper-Fry formula for p t — 2 GeV, Tf = 
.12 GeV , the dash (blue) curve is the Gaussian approximation 
mentioned in the text. 

Now we consider the effect of sound perturbations on 
the particles (pions) spectra. Here we follow the dis- 
cussion of [3] since the same process occurs in our case. 
The perturbations can affect the spectra in two ways. 
We have two variables which have their contributions: 
namely, the flow velocity u M and temperature T. The first 
effect is related to the u p in the exponent of (??) which 
should be corrected by an extra term of the first order due 
to sound. The second effect is related to the first order T 
perturbations. One can consider the so called "zeroth" 
order fireball and the fireball with the "hot spot" induced 
by an initial-state perturbation. In the presence of this 
"hot spot" the corresponding perturbation in tempera- 
ture, ST, is positive and implies a production of extra 
particles as compared to the zeroth order fireball. The 
ST locally delays the freezeout temperature by the for- 
mula 



T (t,r)+ST(t,r)=T f , 



(5.5) 



which also provides an extra volume for containing the 
produced extra particles. Thus the freezeout surface is 
deformed in terms of the temperature and volume. Gen- 
erally, it provides a little larger flow as compared to that 
of the zeroth order fireball. 
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The velocity perturbation, Su^, added to the velocity 
u M gives 



Up -> Uf_i + Su^ 



(5.6) 



In order to find ST and Su we make use of well known 
hydrodynamical relations: 



9 rr.3 

n = ^ T ' 



e = 3p = 'SnT . 



(5.7) 



(5.8) 



The g' is the total degeneracy factor which counts the 
total number of degrees of freedom, summed over the 
spins, flavors, charge (particle-antiparticle) and colors of 
particles. For the Q GP, g' 3 at freezeout. From the 
equations (5.7) and (5.8) one obtains 



ST 



ir 2 h 3 
12s' J T 3 ' 



(5.9) 



In the denominator we can take T 3 

The next step is to use the first order velocity pertur- 
bation from Ref. [9] which is of the form of 



Su 



5e 



(5.10) 



With the corresponding substitutions from (5.7), (5.8) 
and (5.9), the equation (5.10) results in 



Su^ 



T 3 

-Zcs-^STn^ 
1 o 

„ ST ^ 
-oc s — — n? , 



(5.11) 



where nj: is the unit vector in the direction of the mo- 
mentum of phonons which are locared in the points of 
the "crosses" in the transverse xy plane, such as in 



Fig. 2. Thus we have the equations (5.9) and (5.11) which 



determine the first order perturbations of the tempera- 
ture and velocity. The only unknown is the perturbation 
in energy density, 5e, which should be found as well. In 
the appendix we work out the procedure by which the Se 
is determined from the wave equation (see Appendix B 
for more details). 

However, there is also another effect related to the Su 
which is of the first order of perturbation due to expan- 
sion of the fireball: Su = HSr, where r is the radial 
distance from the center of the fireball. From our hydro 
parameterizations we have found an approximate func- 
tion for the temperature (at freezeout time tf = 12 fm) 
as a function of r.It is approximated as 



T(t f ,r) 



-1.75 r 2 

fm 2 



19.20r 
fm 



+ 99.67 jMeV. (5.12) 



It is relevant to note that the Sr gives the increase of the 
fireball's volume (also giving rise to ST). It is equal to 



Sr = 



ST(t,r) 



[dT(t f ,r)/dr]\ rf 



(5.13) 



which, e.g., for 
- 5MeV depend- 



The ST is determined from (5.9) 
dE/dx = lGeV turns out to be 4 -4 
ing on the angle (j>. Then due to the effect of the volume 
change we will have 



Su^ = 



HSr ri? = 

5T(t,r) 



= H \nrru y \ 11 "r, (5-14) 
[dT(tf,r)/dr]\r f 

where np is the unit vector in the radial direction of the 
phonons at the points of the "crosses" in the transverse 
xy plane, such as in Fig. 2. Thus taking into account 
all the above-mentioned effects, for the exponent in the 
Cooper-Frye formula one can write 



T f 



p U pu 



With the perturbations one obtains 



T f 



pu 



pu vSu^ pSvffl 

h h 

Tf T f T f 



(5.15) 



(5.16) 



= Su^ determ ined from (5.11), and 
determined from (5.14). 



where l^u^ 1 ^ 

|^( 2 ) | =SU<- 2 

In heavy ion collisions some amount of secondaries 
originate from jets produced via hard scattering. When 
the transverse momentum p t is large enough, tha hard 
component of it (in the spectrum of particles) "sup- 
presses" the hydrodinamical spectrum, though the jet 
quenching can be significant. So the hydrodynamics does 
not work at large pt's. Such a case occurs also at rela- 
tively large viscous corrections to the viscous term in the 
stress tensor. These corrections get larger with increase 
of pt, and will not be small at some point when one com- 
pares them with the ideal term of the stress tensor. 



VI. COMPARISON TO THE RHIC DATA 

Finally, we can carry out calculations of the pions' 
spectra including the sound wave perturbations, and do 
a comparison with the STAR's two-pion correlation func- 
tions. 

The effect of both components of the flow we called 
above Su^\ Su^ are calculated, for each a, and then 
summed over. We remind that the former effect is the 
effect of perturbed flow, normal to Mach cone, and the 
latter is the effect of extra matter added by the wave, 
which flows together with the rest of matter in the lo- 
cal radial flow direction. The results, for p t = 1 and 
2 GeV/c are shown in Fig. 10 As shown there, the av- 



erage over half plane produce a signal shifted in angle, 
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FIG. 10: (Color online) The spectrum vs <j> p at p t = 1 GeV 
and 2GeF, for dE/dx = 1 GeV/fm. The contribution to 
the spectrum mostly comes from a phonon at the "cross" in 
Fig. 2. The red and blue dashed lines show contributions of 
the jets in the upper (a > 0) and lower (a < 0) half plane in 
Fig.l. 



but in average together with another half-plane it cre- 
ates a plateau-like sum. As p t grows, it starts develop a 
double-hump structure reminiscent of the original Mach 
cone predictions (but with a different angle). 

An example of hard-soft correlation function exhibited 
in Fig.fUl is from STAR publication [33]. The trigger 
Pt = 6.. 10 GeV, which we consider to be well in the hard 
hadron domain. The open (black) circles correspond to 
dAu collision and are shown for comparison only, as in 
this case there is no matter for jets to interact with: they 
display narrow peak at <f> — tt. The central AuAu data 
are shown by (blue) filled circles and open (red) squares, 
for two rapidity slices. Two pictures are for two bins of 
the associate particle (see the captions). 

The overall shape and the width of the angular distri- 



6.0 <p tria < 10.0 GeV/c 




FIG. 11: (Color online) From |33| . Background-subtracted 
azimuthal angle difference distributions for near-central col- 
lisions (fraction of the total cross section 0-12%). The asso- 
ciated particles have the range of pr between 0.5 — lGeV/c 
(upper figure), between 1.5 — 2.5 GeV/c (lower figure), and 
the trigger particles have pr ranging from 6.0 to 10.0 GeV jc. 
The data for Au+Au collisions are shown by the solid cir- 
cles and for d+Au by the open circles. The rapidity range 
is 1 77 1 < 1 and as a result the rapidity difference is |A?j| < 2. 
Open red squares show results for a restricted acceptance of 
|A?7| < 0.7. The solid and dashed histograms show the upper 
and lower range of the systematic uncertainty due to the vi 
modulation subtracted. 



bution seem to be reproduced by the theory. 

While comparing the two, one should however remem- 
ber that: 

(i) our calculation include only energy deposited into the 
media and the sound waves. Sufficiently peripheral jets 
(with a R3 90°) have punch-through remnants of the jets 
which are not included. This component is increasing 
with associate particle pt, filling the place between the 
two humps in the lower picture. 

(ii) theory calculation is done for pions having the same 
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rapidity as the associate jet, while the data are for cer- 
tain rapidity interval. However, as comparison of the 
solid (blue) circles and open (red) squares shows, the dif- 
ference in shape is not too strong. This issue definitely 
need to be studied more, to separate what we call "cone" 
from the rapidity-independent "circles" . 
(iii) theory calculation have no fluctuations in dE/dx, 
and thus assume completely deterministic scenario in 
which all jets stop at some fixed points, when its energy 
is all lost. Perhaps in reality fluctuations are not small, 
and there are punch-through jets at any a. 

Completing the paper we emphasize that new data 
from LHC have well reconstructed di-jet events, in which 
transverse momenta and rapidities of both trigger and 
associate jets are detected. Soft particles angular dis- 
tributions, both in azimuth and rapidity, will obviously 
provide much better understanding of the issues involved. 
However current statistics of such events ~ 1000 is per- 
haps too small to address it. We will report calculations 
of the sound wave contributions in such conditions else- 
where. 
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VII. APPENDIX A: THE FREEZEOUT 
SURFACE 

We approximate the freezeout surface as the isoterm, 
with Tf = 120 MeV. The parameterization of it is based 
on the results of other authors, namely Refs. [T7]. [25] 
and |26j . The surface we use is parameterized as 



r(r) 



-0.0011 * r 5 .995 + 0.0266 * r 4 . 99762 (7.1) 



-0.2367 * r 4 + 0.9772 * r J 
+0.7* r + 17.657 



1.9095 *r z 



In [17], it was noted that the freezeout time changes 
only a little if one goes from SPS (tf — 10 fm) to 
RHIC {t f = 11 fm). In [23], the freezeout time 
for central collisions is calculated to be « 12 fm (at 
Tf = 130 fm forRHIC)In [26], the value obtained for the 
freezeout time is tf = 12 fm/c, however, the tempera- 
ture distributions can have sizable spreads in the range 
of - 100-160MeF. In [27], having the results of [25]. it 
was pointed out that the evolution time to kinetic freeze- 
out - say until T f w UOMeV is tf ~ 12 - Ufm/c in 
both cases of longitudinal and transverse dynamics; the 
timescale which is most directly probed by femtoscopy. 
Nonetheless, from all these scenarios, having a worked 
out combined scenario from [T7] and [23] at our dis- 
posal, we take the aforementioned values for the freeze- 
out time and "kinetic" temperature as averaged values 
for our ansatz: namely < tf >= tf — 12 fm and 



< Tf >= T 



f 



120 fm. Such a combined scenario is 



Usually one separates "kinetic" and "chemical" freeze- 
outs in which elastic and inelastic scattering rates are 
involved. Since different secondaries (pions, K mesons, 
nucleons, J/tp particles, etc) in fact have quite differ- 
ent elastic cross sections, the "kinetic" surfaces should 
be different for each species. We do not discuss all such 
complications in our work, and consider only one type of 
secondaries, the pions, the spectra of which are shown in 
the last section. 



VIII. APPENDIX B. NORMALIZATION OF 
THE MACH CONE WAVE 

Suppose that the relative changes of the density p and 
pressure p in the medium of the fireball are much small 
as compared to their unperturbed values: Sp -C po and 
Sp -C pa such that 



p = po + Sp and p = p + Sp . 



(8.1) 



Suppose that this is also the case for the velocity of par- 
ticles of the medium: u <C c s . 

Then having the continuity equation and Euler's equa- 
tion 



dp 
dt 
du 
dt 



+ Vpu = 
+ (mV)« 







p 



(8.2) 



one can rewrite this system as 

d(S P ) 
dt 

du V{Sp) 
dt po 



Po V • u = 



= 0. 



If now we proceed to relativistic notations Sp 
Po = e o+Po> then the system i 
form: 



(8.3) 



Se and 



3 1 will have the following 



d(Se) 
dt 

{to+Po) 



+ (eo +Po)V • u 

du 
dt 







V(Sp) = 0. 



(8.4) 



Actually, here the first equation is the energy conserva- 
tion and the second one is the Newton's second law. 

The small change in pressure p is related to the small 
change in energy density e by the adiabatic equation: 



Sp 



dp 



Sp = czSe. 



(8.5) 



The standard procedure with the system (8.4) leads to 
the wave equation both for Sp(x, y, z, t) and Se(x, y, z, t) 
(see Ref.): 



based on the joint consideartion of these papers. 



V 2 (Sp) 



1 d 2 (S P ) 
c? dt 2 



(8.6) 
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V 2 (<5e) 



1 3 2 {5e) 
c! dt 2 



0. 



(8.7) 



We focus on the last equation which is a second or- 
der partial differential homogeneous equation. In case of 
a jet propagation which is a source of sound waves, in 
the rhs of (8.7) there will be some function f(x,y,z,t) 
describing the source, and making the equation nonho- 
mogeneous: 



V 2 (^) - 



1 d 2 {8e) 
£ dt 2 



f(x,y,z,t) . 



(8.8) 



Our goal is to solve the equation (8.8 1 which is carried 



out as follows. First we should find the solution of the 
homogeneous equation (8.7). In turn it can be rewritten 



in terms of cylindrical-coordinate wave equation such as 



1 8 2 (5e) _ 

c 2 dt 2 

d 2 {6e) 1 d(6e) 1 d 2 (8e) d 2 (5e) 



dr 2 



r t dr t 



dip 2 



dx 2 



, (8.9) 



where Se = 5e(rt,ip,x,t), r\ 



V 



z 2 , and the ip is 



the azimuthal angle. The solution to the cylindrical- 
coordinate wave equation can be found using the sep- 
aration of variables: 

Se(r t , ip, x, t) = R{r t ) <f>(ip) X(x) T(t) . (8.10) 

Solving by this way one gets two final solutions in what 
follows. 

6e ki ,k, n = A± k n J n (k t r t ) e ±*v^F *±*c.t±*w , (8 . n) 
<W,« - B^ n Y n {k t r t ) e ±iV^^c 3 t±in V 5 (8 J2) 



where k t and k x = \Jk 2 — k 2 are the transverse and 
longitudinal momenta, and J n [ktTt) & n d Y n (k t r t ) are the 
Bessel functions of the first and second kind. The k n 

and B kt k n are some constants coming from solutions 
with the separate variables. The notation in the ex- 
ponent of ( pTTT] ) (or ( |8l2| ) means that there are eight 
combinations of the sum corresponding to positive and 
negative values of x, k and ip. At azimuthal symmetry 



that n = in (8.11) and" 
them as 



the term ^ d from (8.9 ) is dropped out which means 



(8.12). Thus we can rewrite 



<fe*i,* = A i,kMhn) e ±<V^l*+^) , (8.13) 
fek t ,k = Bt k Y (k t r t ) e ±HV^Ix+^ , (8.14) 



In this case in the exponent of (8.13 1 (or (8.14)) there are 



four combinations of the sum corresponding to positive 
and negative values of x and k. 



We solve this last equation using (8.13) and/or (8.14). 
Suppose we use (8.13). Then the operators in the l.h.s. 



of (8.15) take the following forms: 



d 2 (Se) 
dt 2 



A^ k Uk t r t )e ± ^V^I-+^) , (8.16) 



-c 2 y 2 (Se) = 
-- c 2 s (k 2 + k 2 x ) x 

Ai. k Jo(k t r t )e ±l (V^!*+^) 



117) 



By means of these two expressions one can rewrite the 
solution of (8.15) in the u and k space, such as 



[-lo 2 + c 2 s {k 2 t + hi)] <j> uk = f uk = f(u, k) . (8.18) 

The /(w, k) is determined by the Fourier transform of 
f(x, y, z, t) which one can choose as a function of the 
form (l/ct f )5(±x - ct)5(y)S(z): 

f(u,k) = 

dtdxdydze ±l ^ k2 ~ k t x+ult) x 

1 



x — S(±x - ct)S(y)S(z) 

Ctf 



(8.19) 

The c is the speed of light. From (8.18) the (j> uk is repre- 



sented as 



z 2 k 2 



(8.20) 



Plugging the result of the integral of ( 8.19 ) into this equa- 
tion we obtain 



>uk)+ 



2tt 



c 2 k 2 



61" 



2tt 



- (-M) • (8-21) 



The (<^jfc)+ is for positive k x , and the ((j) u k)— is for 
negative k x . Carring out the inverse Fourier transform 
of <j>uk along with ( 8.21| ) we get 

[<f>(t,x,y,z)] + = [5e(t,x,y,z)} + = 



and (8.8) as 

d 2 (8e) 

dt 2 



c 2 V 2 (Se) = f(t,x iy ,z). 



3.15) 



[<j>(t,x,y,z)]- = [Se(t,x,y,z)]- = 

d ^0e^-^\^ k )^ (8.22) 
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which consecutively results in 



Hereby, we further have 



[Se{t,x,y,z)] + 
1 



^i(k x ct— (k x x-\-k y y-\-k z z)) 



cm+k 2 y + kf) - C 2 k 2 > 



dkr^ dkr^ 



[5e{t,x,y,z)] 
{2m) 



{2-Kfctf 

{kt/cMhn) 



2k t 



,j 2 -e 2 

«t=±i«/ — k x 



then 



[Se{t,x,y,z)]- 
1 



e Ti(~k x ct+(k x x+k y y+k z z)) x 



{2-Kfctf 

dk 

r 2(U2 i U2 I 1.2 \ _ r 2h.2 ' 

c s \ K x ' K y ' K z) c K x 



dfcj£ dky dkr^ 



(8.23) 



[6e{t,x,y,z)] + = 
1 



{2n) 3 ct f . 
x / dk t dtp 



jU e Tik t r t cos ip 



.2 / / 2 _i_ ^.2j 



S 2 fc 2 ' 



[<5e(i,x,y, z)]_ = 
= (2ttz) /" eT4 (-fe, ct+ fe xa; ) 
{2w) 2 ct f J x 
-{k t /c 2 s )Mk t r t y 



2k t 



[5e(t,x,y,z)} + = 



, (8-27) 



4irc 2 .ct 



dk T eT^ct-k^x) x 



x J ( ±M/ Cs n° nk x ) , 



[Je(t,a;,y,z)]_ = 

— \ I Ah pTii-kxCt+k^x) 

-j2n)^r f J dkxe 

/fc e ±ik t r t cos ip 
dhdip c 2 M+w-*kr (8 - 24) 



[Se{t,x,y,z)\. 



4irc 2 ctf 
x J ±i\ 



dk x e Tl( -- k * ct+k * x) x 



(8.28) 



[Se{t,x,y, z)] + 
1 



(2n) 2 ct f 
x / dk t 



dk x e Ti(k * ct - k * x) x 



(fc t /c*) Jo(fctrt) 



k 2 + [{c 2 -c 2 )/cl]kr 



[Se{t,x,y,z)]- = 
1 



(27r)2 rf/ 



dk efii-kxCt+k^x) x 



x / dkt 



(k t /c 2 ) J Q {k t r t ) 
k 2 + [(c 2 ~c 2 )/c 2 ]k 2 x 



In the integral over fc t the first order pole is 



k, = ±it 



(8.25) 



(8.26) 



Ultimately, the integration of two integrals in (8.28) 



gives the following solutions for the function 5e(t, x, y, z) 

[6e(t,x,y,z)} + = 
= 0, 

1 1 



and 



2irc s ct f ^Jc 2 {x - ct) 2 + (c 2 - c 2 )r 2 ' 
1 1 
2nc s ct f ^/cf( x - ct) 2 + (c 2 - c 2 s )r 2 ' 
0, (8.29) 



[6e{t,x,y,z)}_ = 
= 0, 

1 



1 



2nc s ct f y / cf(x - ct) 2 + (c 2 - c 2 )rf 

1 1 
~ 2nc s ct f ~f<? a { x - ct) 2 + (c 2 - c 2 )r 2 ' 
0. (8.30) 
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Combining the solutions which are in (8.29) and (8.30) 
we obtain 

[8e(t,x,y,z)} + = [Se(t, x, y, z)]- = 

= ±^ I . 

2irc s ct f ^/c 2 (x - ct) 2 + (c 2 - c 2 )(y 2 + z 2 ) 

= 0. (8.31) 



Generally, the function f(t,x,y,z) in equation (8.15) 



can include some constant A which has a dimension of 
GeV/fm. Taking also into account this constant we have 



5e(t,x,y,z) 
1 



{8.32) 



2TTC s ct f ^c^Jx - ct)2 + ( c 2 _ c 2)( y 2 + z 2-j 

One can make this constant A equal to dE/dx which is 
the loss of the jet energy per unit path length. 
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